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I. 


SUMMARY 


In  addition  to  the  previously  developed  partially  cavitati  g 
cascade  theory,  two  new  flow  models  were  constructed  in  search 
of  a  better  flow  model  for  determining  accurate  force  coef¬ 
ficients.  Effort  has  been  made  for  obtaining  i)  physically 
acceptable  flows,  particularly  the  location  of  cavity  boundary 
and  ii)  smooth  matching  of  the  flow  characteristics  between  the 
partially  cavitating  and  supercavitating  flow  regimes.  Based 
on  the  numerical  results  made  with  these  flow  models  for 
practical  blade  profiles  taken  after  a  supercavitating  pro¬ 
peller  it  was  found  that  no  single  flow  model  developed  above 
could  handle  the  complete  set  of  cascade  geometries  and  inci¬ 
dence  angles.  One  theory  was  supplemental  to  the  other  and  no 
definite  guideline  was  discovered  for  selection  of  an  appropriate 
flow  model  for  a  specified  flow  condition  to  be  solved  except 
for  a  few  weak  evidences.  The  data  calculated  with  one  of  the 
above  theories  were  used  in  the  high  speed  propeller  theory 
for  predicting  the  performance  of  a  supercavitating  propeller 
operating  at  the  partially  cavitating  regime.  The  results 
showed  a  slight  improvement  over  the  previous  ones,  indicating 
a  good  possibility  of  increasing  the  prediction  capability  of 
the  propeller  theory  for  the  partial  cavitating  flow  regime. 

In  doing  this,  a  more  accurate  theory  will  have  to  be  developed. 
Before  making  further  effort  on  developing  such  a  theory,  many 
unknown  aspects  regarding  the  boundary  condition  of  the  par¬ 
tially  cavitating  cascade  flow  should  be  clarified.  A  systematic 
experimental  study  for  such  flow  is  now  in  order  and  strongly 
recommended . 
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1. 


INTRODUCTION 


A  performance  prediction  method  for  cavitating  marine  propellers 
has  been  developed  by  combining  the  propeller  lifting-line 
theory  with  the  supercavitating  cascade  theory  (Furuya  1976, 
1978b) .  This  method  was  applied  for  calculating  the  performance 
characteristics  of  a  few  supercavitating  propellers  such  as 
TMB  Model  3770  and  Hydronautics  7607.02.  It  was  found  that 
the  predicted  force  coefficients  with  the  method  compared 
favorably  with  their  experimental  data  over  a  range  of  the 
advance  speed  J  for  which  the  supercavitating  condition  prevails. 
On  the  other  hand,  for  larger  J's  for  which  the  partially  cavi¬ 
tating  flow  appears  on  a  part  of  or  whole  propeller  blade,  the 
predicted  thrust  and  torque  coefficients  rapidly  departed  from 
those  of  experiments  (see  Figures  2.12  and  2.14,  reproduced 
from  Furuya  1978b) .  It  was  considered  that  this  discrepancy 
was  attributable  to  the  use  of  the  erroneous  supercavitating 
cascade  data  in  the  propeller  theory  for  the  partially  cavi¬ 
tating  flow  regime. 

The  objective  of  the  whole  project  was  therefore  first  to 
develop  the  partially  cavitating  cascade  theory  and  then  to 
incorporate  the  data  calculated  with  such  theory  into  high¬ 
speed  propeller  theory.  During  the  initial  phase  under  the 
GHR  project  the  backbone  for  the  partially  cavitating  cascade 
theory  was  established  (Furuya  1978).  Computer  programs  (PCAS 
and  PCASL)  developed  based  on  the  above  mathematical  formula¬ 
tion  were  applied  to  a  cascade  of  plano-convex  blades  with 
various  blade  thickness.  Representative  numerical  calcula¬ 
tions  made  with  these  programs  showed  a  good  agreement  with 
experimental  data  of  Wade  and  Acosta  (1967) . 

The  remaining  tasks  for  the  present  phase  were  to  generate 
sectional  data  for  practical  blade  profiles  of  a  supercavi¬ 
tating  propeller  with  the  previously  developed  cascade  theory 


1 


and  then  to  use  such  data  in  the  propeller  program  for  pre¬ 
dicting  the  propeller  performance  operating  under  the  partially 
cavitating  regime.  When  the  blade  sectional  profile  of  Hydro- 
nautics  7607.02  supercavitating  propeller,  having  rather  large 
camber,  was  used  in  PCAS  and  PCASL,  a  problem  of  the  calculated 
upper  cavity  boundary  crossing  the  blade  surface,  which  is 
physically  not  acceptable,  arose.  This  was  not  surprising  since 
the  flow  model  used  here  had  no  control  on  the  location  of  the 
cavity  closure  point  and  this  type  of  problem  is  quite  commonly 
experienced  in  supercavitating  blade  design  procedure  (see, 
e.g.  Yim  and  Higgins  (1975)).  In  order  to  overcome  this  prob¬ 
lem  we  developed  a  new  flow  model  for  the  partially  cavitating 
cascade  flow  in  which  the  boundary  condition  for  the  cavity 
end  point  to  stay  on  the  upper  blade  surface  was  added  but 
with  allowing  the  wake  to  open  at  downstream  infinity.  In  the 
following  section  the  theory  formulation  for  this  new  flow 
model  and  computation  for  the  sectional  blade  profiles  of 
Hydronautics  7607.02  propeller  as  well  as  the  application  of 
the  data  obtained  there  for  the  propeller  theory  will  be 
described  in  detail. 

As  will  be  seen  in  Section  2,  this  new  flow  model  was  found  to 
be  not  quite  satisfactory  in  accuracy  for  the  data  obtained. 

The  lift  and  drag  forces  calculated  from  this  partially  cavi¬ 
tating  cascade  theory  sometimes  matched  smoothly  with  those  of 
the  supercavitating  cascade  theory  as  £c/c  approached  unity, 
but  sometimes  not,  depending  upon  the  cascade  geometry  and 
incidence  angle.  It  was  discovered  that  this  mismatching  was 
particularly  eminent  for  rather  small  incidence  angle  cases. 
Close  investigations  on  the  flow  models  used  for  both  super¬ 
cavitating  (s/c)  and  partially  cavitating  (p/c)  theories  as 
well  as  the  calculated  results  revealed  two  possible  reasons 
for  the  discrepancy.  First,  the  base  pressure  between  the  s/c 
and  p/c  theories  is  totally  different  no  matter  how  closely  the 
length  of  cavity  approaches  to  that  of  chord;  it  is  pc  for  the 
former  whereas  it  is  p2  for  the  latter.  Secondly,  it  was 


shown  in  the  results  obtained  with  the  s/c  theory  that  the 
cavity  thickness  was  so  thin  that  it  either  barely  cleared  the 
upper  blade  surface  or  crossed  it  for  small  incidence  angle 
cases.  There  exists  no  remedy  for  the  cavity  boundary  of  the 
s/c  flow  crossing  the  blade  surface.  It  simply  means  that  the 
data  for  such  cases  must  be  ignored  since  such  flows  should 
not  physically  exist. 

For  a  further  improvement  on  the  partially  cavitating  cascade 
theory  the  third  flow  model,  i.e.,  "double  wake"  model,  was 
introduced  here.  In  this  flow  model,  the  near  wake  and  far 
wake,  one  after  another  immediately  behind  the  blade  having 
blunt  trailing  edge,  were  used  in  the  same  way  as  is  observed 
in  the  actual  wake  flow.  In  choosing  the  pressure  for  the 
near  wake.  Professor  Acosta  at  the  California  Institute  of 
Technology  suggested  a  linear  interpolation  between  those 
inside  the  cavity  and  far  wake.  Although  the  validity  of  this 
interpolation  method  should  await  comparison  with  experimental 
data,  the  method  at  least  guarantees  the  smooth  blending  of 
the  base  pressure  between  the  two  flow  regimes.  In  Section  3 
this  double  wake  model  will  be  described,  followed  by  repre¬ 
sentative  numerical  results. 
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2.  MODIFICATION  OF  PARTIALLY  CAVITATING  (p/c)  CASCADE 

THEORY  FOR  CAMBERED  FOILS 

2 . 1  THEORY 

The  nonlinear  theory  for  partially  cavitating  cascade  flows 
previously  developed  (Furuya  1978)  was  applied  to  the  blade 
profiles  of  Hydronautics '  7607.02  propeller.  Although  convergent 
solutions  were  obtained  with  this  theory,  we  encountered  a  prob¬ 
lem  of  the  cavity  end-point  crossing  the  foil  upper  surface  and 
ending  at  a  physically  unacceptable  point.  Figure  2.1  shows 
one  such  example.  The  blade  profile  and  cascade  geometry 
used  for  this  calculation  were  those  at  50%  radial  station  of 
7607.02,  incidence  angle  =  4°  and  %c  =  .35.  This  type  of 
problem  is  well-known  in  the  design  method  of  supercavitating 
hydrofoils  or  cascades  (see  Yim  and  Higgins  1975)  where  the 
pressure  distribution  which  provides  a  good  lift-to-drag 
ratio  is  usually  specified  and  the  corresponding  foil  pro¬ 
file  should  be  determined.  In  many  cases  designers  will  find 
that  the  calculated  upper  cavity  boundary  crosses  over  the 
blade  surface. 

In  order  to  overcome  this  problem,  a  modification  has  been  made 
on  the  existing  theory  (Furuya  1978)  and  then  the  computer  pro¬ 
gram;  an  additional  constraint  has  been  placed  to  force  the 
cavity  end-point  to  close  on  the  upper  surface  of  the  foil. 

Along  with  this  additional  boundary  condition,  a  new  unknown 
variable,  s,  was  introduced  to  accommodate  the  additional  equa¬ 
tion  where  e  is  taken  to  be  the  ratio  of  the  blade  spacing  to 
the  flow  width  at  downstream  infinity.  It  means  that  with  the 
additional  condition  forcing  the  cavity  end-point  to  end  on  the 
body  boundary  it  is  no  longer  possible  for  the  wake  to  close  at 
downstream  infinity  (see  Figure  2.2).  This  extra  condition  and 
unknown  parameter  therefore  will  change  the  flow  model  from  the 
closed  wake  model  to  the  open  wake  model.  It  is  noted  that  our 
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open  wake  model  is  different  from  Wu's  open  wake  model  (1968) 
since  the  pressure  inside  the  wake  in  this  model  is  a  constant 
everywhere  and  equal  to  that  of  the  downstream  infinity.  Again 
two  types  of  the  boundary  value  problems,  one  specifying  the 
cavitation  number  (computer  code  "PCASE")  and  the  other  specify¬ 
ing  the  length  of  cavity  (computer  code  "PCASLE"),  were  formu¬ 
lated  and  the  computer  codes  were  tested  in  several  cases. 
Although  most  of  the  boundary  conditions  remain  the  same  as 
those  of  the  previous  theory  (Furuya  1978) ,  the  new  boundary 
value  problem  used  here  will  be  described  in  the  complete  form 
for  a  clear  understanding  of  the  theory  change. 


The  solution  form  for  the  hodograph  variable  w  also  remains  the 
same  (see  Furuya  (19  78)), 


The  unknown  parameters  in  the  present  problem  include  U2 ,  a2 
A,  b,  c,  f  and  e,  requiring  seven  equations  for  the  unique 
determination  of  these  parameters.  The  boundary  conditions  to 
be  used  here  are  described  as  follows: 

(i)  At  the  upstream  infinity, 

oj  ( +  i  £n(l/U2):  2  equations 

(ii)  At  the  downstream  infinity, 

(°°)  *  a2  5  1  equation 

(iii)  Length  of  first  arc  =  S^:  1  equation 

(iv)  End  point  of  the  streamline 
(£»f)  matches  the  physical 
trailing  edge:  1  equation 


(v)  Continuity  equation: 


1  equation 


where 

sgU>  -  j  1  '  «  *  0 

1-1  ,  5  <  0  . 

The  arc  length  condition  is  therefore  satisfied  by  the  follow¬ 
ing  equation; 

f 4=  S]_  -  s1(-l)  =  0  (2.9i 

where  dW/d£  is  given  by  Equation  (4) . 

iv)  End  Point  of  the  Streamline  at  •  =  f  Matching  the  Trailing  Edge 
of  the  Blade: 

Similar  to  the  above  derivation,  the  equation  f_  is  obtained  as 
follows : 


f5  -  S2  “  s  2  (f)  =  0 


(2.10) 


where 


f  r  \ 


In  order  to  determine  the  second  arc  length  S2,  the  end  point  of 
cavity  should  be  known.  The  cavity  shape  is  expressed  in  terms  of 
5 -coordinate? 


XB'  yB  are  the  Physical  coordinates  of  the  blade  leading  edge. 

The  end  point  of  the  cavity  (x  ,  y  )  is  obtained  by  setting  the 

c  c 

upper  bound  of  the  above  integrals  (2.13)  and  (2.14)  to  be  c. 
Once  (x  ,y  )  are  determined,  the  second  arc  length  S,  is  readily 
calculated  by  following  the  body  profile. 


The  continuity  equation  applied  here  at  upstream  and  downstream 
infinities  becomes 


fg  =  cos  (oi1+Y)  •d1  -  U2  cos(a2+Y)*d2  -  0  (2.1 

where  £  is  defined 


vi)  Cavity  Closure  Condition 

The  new  boundary  condition  to  be  added  is  written 

f 7  -=  ^c  *  W  "  0  (2-J 

c 

where  x  =  ^  f  cos  gc(^')  cU '  (2.] 

b 


C 

yc  =  yb  +  ==  /  sin  gcu,)  snr  dr  (2‘2 


f  =  the  upper  surface  shape  of  the  foil. 
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NUMERICAL  RESULTS  WITH  THE  NEW  THEORY 


The  two  computer  programs  (PCASE  and  PCASLE)  developed  above 
were  used  to  calculate  the  lift  and  drag  forces  for  the  propeller 
blade  sectional  profiles  of  the  Hydronautics  7607.02  propeller. 

The  previous  problem  of  the  cavity  end-point  crossing  the  blade 
boundary  has,  of  course,  been  resolved  with  these  new  programs. 
Figure  2.3  shows  the  calculated  cavity  boundary  for  exactly  the 
same  case  as  that  shown  in  Figure  2.1,  proving  that  the  new 
boundary  condition  for  the  cavity  end  point  has  properly  been 
implemented  in  the  theory  and  computer  program. 

Although  the  cavity  boundary  problem  was  solved,  a  new  type  of 
problem  has  emerged.  Figures  2.4  to  2.9  show  the  lift,  drag 
coefficients  and  cavity  length  as  a  function  of  cavitation  num¬ 
ber  over  supercavitating  and  partially  cavitating  regimes. 

The  force  curves  of  solid  lines  in  these  figures  are  those 
which  had  been  calculated  with  the  supercavitating  (s/c)  cascade 
theory  (Furuya  1975)  ,  the  same  as  those  shown  in  the  previous 
report  {Furuya  1978b) ,  whereas  those  of  broken  lines  are  the 
newly  calculated  values  with  PCASE  and  PCASLE,  covering  the 
partially  cavitating  regime.  It  is  seen  from  Figures  2.4  and 
2.7  that  the  lift  curves,  as  they  approach  a  point  of  l /c  =  1 
from  the  two  regimes ,  in  some  cases  match  but  in  other  cases 
do  not,  depending  upon  the  incidence  angle  and  spanwise  loca¬ 
tion.  On  the  other  hand,  the  drag  values  calculated  from  the 
partially  cavitating  cascade  theory  did  not  agree  with  those  of 
the  s/c  cascade  theory  near  l  / c  *  1  for  all  cases.  Table  2.1 

G 

shows  most  of  those  data  near  l  /c  =  1. 

c 

Figures  2.6  and  2.9,  which  show  the  relationship  between  the 
cavity  length  and  cavitation  number,  more  clearly  indicate  the 
mismatching  between  the  two  theories.  In  addition  to  the  mis¬ 
matching  of  the  data  at  £  /c  =  1,  we  ran  into  difficulty  in  under- 

C 

standing  the  physical  meanings  of  the  computed  results  for  smaller 
incidence  angles  than  =  4°,  such  as  =  2°,  0°,  etc.  For 
those  smaller  incidence  angles,  the  cavitation  number  which 
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was  determined  in  the  program  decreased  as  the  specified  length 
of  cavity  was  decreased.  It  must  be  pointed  out  that  we  shall 
not  totally  blame  the  partially  cavitating  cascade  theory  for 
the  mismatchings  discovered  above.  When  the  upper  cavity 
boundaries  calculated  by  the  s/c  cascade  theory  were  plotted, 
it  was  found  that  they  crossed  the  blade  boundaries  particularly 
for  small  incidence  angle  cases  (see  Figures  2.10  and  2.11) 
at  the  spanwise  locations  50%  and  60%,  respectively.  In 
Figure  2.10  at  ot^  =  8°  the  calculated  cavity  boundaries  cleared 
the  upper  trailing  edge  of  the  blade  for  all  o's  whereas  at 
*  4°  the  last  o  for  which  the  cavity  boundary  cleared  the 
edge  was  at  .22.  It  means  that  the  supercavitating  cascade 
data  at  =  4°  for  a's  less  than  about  .22  are  physically 
unacceptable  for  this  blade  profile  at  x  =  .5.  A  similar  situ¬ 
ation  also  exists  in  Figure  2.11  for  the  flow  at  the  spanwise 
location  of  60%. 


13 


2.3  APPLICATION  OF  THE  NEW  p/c  CASCADE  DATA  TO  A  CAVITATING 

PROPELLER 

The  objective  of  the  present  work  is  to  determine  the  effect 
of  using  the  data  of  the  partially  cavitating  cascade  flow 
on  the  propeller  performance  prediction  method.  Previously, 
even  for  the  regime  of  the  partially  cavitating  condition,  the 
data  of  supercavitating  cascade  theory  were  used,  leading  to 
significant  deviation  of  the  predicted  performance  from  the 
experimental  results  (see  Furuya  1978b) . 

Although  the  discrepancies  existed  between  the  s/c  cascade 
data  and  the  present  p/c  cascade  data,  these  partially  cavi¬ 
tating  cascade  data  were  read  from  Figures  2.4,  2.5,  2.7,  and 
2.8  and  fed  into  the  propeller  program  "SCPR0P2"  of  Furuya 
(1978b)  as  part  of  the  input  data.  We  know  in  advance,  there¬ 
fore,  that  the  calculated  results  from  SCPR0P2  might  not  be  as 
accurate  as  they  should  be  but  would  indicate  the  validity  of 
the  present  concept  of  incorporating  the  partially  cavitating 
cascade  theory  into  the  cavitating  propeller  theory. 

Calculations  were  made  for  two  different  cavitation  numbers, 
aVa  =  .343  and  3.0  for  the  Hydronautics  7607.02  propeller  (°Va 
is  the  cavitation  number  based  on  the  ship  advance  speed) . 

The  results  for  the  thrust,  torque  coefficients  and  efficiency 
are  shown  in  Figures  2.12  and  2.13,  respectively,  in  comparison 
with  the  experimental  data  of  Peck  (1977)  and  Bohn  (1977)  as 
well  as  with  the  previous  theoretical  data  of  Furuya  (1978b) . 

For  aVa  =  .343  which  is  the  design  cavitation  number,  the  new 
calculations  provided  a  slight  improvement  over  the  previous 
prediction,  but  not  significantly.  In  Figures  2-13 (a-d)  the 
.detailed  information  for  the  local  incidence  angles,  cavitation 
numbers  and  lift  and  drag  coefficients  is  shown  as  a  function 
of  spanwise  location  x  at  various  J's.  It  is  seen  from  these 
data  that  at  J  =  0.9  to  1.2,  the  local  flow  remains  in  the  super¬ 
cavitating  condition.  The  effect  of  using  the  partially  cavi- 
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tating  cascade  data  in  the  propeller  program  is  not  clearly 
seen  except  that  the  calculations  of  force  coefficients  were 
influenced  by  the  lower  values  of  p/c  cascade  data  during  the 
force  interpolation  procedure.  Furthermore,  it  should  be  noted 
that  for  the  cases  of  small  incidence  angles  the  p/c  cascade 
theory  had  difficulty  in  obtaining  convergent  solutions.  Since 
such  data  were  scarce,  "SCPR0P2"  extrapolated  from  those  of 
higher  incidence  angles.  This  may  be  one  of  the  reasons  for 
which  the  performance  prediction  of  the  propeller  at  partially 
cavi tating  regime  with  the  present  method  has  not  provided 
significant  improvement. 

On  the  other  hand,  for  oVa  =  3.00,  the  present  method  has  over¬ 
corrected  the  thrust  and  torque  coefficients  for  a  range  of 
J's  in  which  the  partially  cavitating  condition  appears  (i.e., 

J  .55),  see  Figure  2.14.  From  the  detailed  information  in 
Figure  2.15  a  sudden  jump  in  various  key  quantities  such  as 
a&,  CL,  etc.,  from  J  =  . 50  to  J  =  . 60  is  seen,  indicating  that 
"SCPROP2 "  program  used  the  data  of  partially  cavitating  range 
for  J's  larger  than  0.60.  An  abrupt  change  in  and  in 
Figure  2.14  apparently  shows  use  of  inappropriate  data  for  the 
current  propeller  analysis.  As  has  been  mentioned  at  the 
beginning  of  this  section,  it  seems  that  the  data  obtained 
with  the  modified  partially  cavitating  cascade  theory  are  not 
accurate  enough  to  be  used  for  the  propeller  performance  pre¬ 
diction  at  the  partially  cavitating  regime. 

It  can  be  readily  concluded  from  the  above  findings  that  more 
accurate  data  for  the  partially  cavitating  cascade  flow  are 
necessary.  It  is  for  this  reason  that  the  double  wake  flow 
model  has  been  introduced  for  the  partially  cavitating  cascade 
theory  and  this  will  be  described  in  the  following  section. 
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3. 


DOUBLE  WAKE  MODEL  FOR  p/c  CASCADE  FLOWS 


3 . 1  THEORY 

As  has  been  discussed  in  the  previous  section,  there  has 
existed  the  discrepancy  in  forces  and  cavitation  numbers  between 
the  supercavitating  cascade  theory  (Furuya,  1975)  and  present 
partially  cavitating  cascade  theory  as  the  cavity  length  approaches 
that  of  the  chord.  This  posed  a  problem,  particularly  for 
foils  having  blunt  trailing  edges  but  not  for  the  sample  case 
shown  in  Furuya  (1978)  in  which  the  plano-convex  cascade  having 
a  sharp  trailing  edge  was  used.  Figure  3.1  shows  two  different 
flow  configurations  as  the  cavity  length  approaches  the  chord 
length  both  from  the  partially  cavitating  condition  side  and 
supercavitating  condition  side.  The  major  difference  between 
these  two  limiting  cases  lies  in  the  base  pressure;  it  is  p2 
for  the  partially  cavitating  flow  model  whereas  that  is  pc 
for  the  supercavitating  flow  model.  These  two  pressures  are 
quite  different,  pc  being  usually  much  lower  than  p2 .  This 
discrepancy  will  remain  in  these  flow  models  indefinitely  no 
matter  now  closely  the  cavity  length  is  forced  to  become  that 
of  chord.  The  pressure  distributions  shown  in  Figures  3.2(a) 
and  3.2(b)  clearly  indicate  the  above  discrepancy;  particularly 
note  the  differences  in  the  base  pressure. 

Another  aspect  to  be  questioned  regarding  the  flow  model  cur¬ 
rently  used  in  the  partially  cavitating  cascade  flow  is  the 
single  wake  flow  configuration.  In  the  fully  established  wake 
flow  around  a  body  having  blunt  trailing  edge,  two  different 
types  of  wake  regions  are  observed  experimentally.  The  first 
region  appears  right  behind  the  body  and  is  enclosed  by  the 
thin  and  rather  steady  free  shear  layers ,  which  may  be  called 
"the  near  wake".  The  pressure  inside  this'  near  wake  is  nearly 
constant,  designated  pw»  Further  downstream,  the  shear  layer 
gradually  becomes  wider  as  the  vorticity  diffuses,  generating 
an  irregular  flow  and  then  turbulent  mixing  flow.  The  flow 
inside  this  region  is  unsteady  and  stationary  only  in  the  mean. 
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This  portion  of  the  wake  may  be  called  "the  far  wake",  along 
which  the  mean  pressure  gradually  recovers  from  pw  to  p2  of 
downstream  infinity  (p2  *  p^  for  the  unbounded  flow) . 

In  light  of  the  observations  on  the  real  wake  flow  discussed 

above,  a  new  flow  model  for  the  partially  cavitating  cascade 

problem  is  proposed  and  discussed  herein.  The  idea,  which  has 

come  out  of  one  of  the  author's  recent  discussions  with 

Professor  Acosta  at  the  California  Institute  of  Technology, 

is  to  represent  the  real  wake  flow  as  closely  as  possible  by 

using  two  wakes,  behind  the  body,  see  Figure  3.3(a).  In  this 

flow  model,  which  we  call  "double  wake"  model,  the  pressures 

p..  in  the  "near  wake"  and  p-  in  the  "far  wake"  are  assumed 
w  z 

constant.  The  problem  of  this  flow  model  is  the  method  of 
determining  p  ,  which  is  expected  to  resolve  the  discrepancy 
which  existed  in  the  previous  flow  model.  Professor  Acosta's 
proposal  goes  to  this  point;  i.e., 


It  means  that  the  "variable"  pressure  pw  in  the  near  wake  is 
determined  depending  upon  the  cavity  length.  As  the  cavity 
length  becomes  that  of  chord,  pw  becomes  the  cavity  pressure 
Pc  whereas  pw  =  p2  is  used  when  the  partial  cavity  totally 
disappears.  The  most  significant  feature  of  this  variable 
pressure  concept  is  that  the  smooth  blending  of  the  partially 
cavitating  flow  into  the  supercavitating  flow  may  be  obtained. 
Although  the  validity  of  linearly  interpolating  p  between  p 

W  v 

and  p2  should  be  proved  in  comparison  with  experimental  data, 
Equation  (3.1)  holds  true  at  least  at  the  limiting  case  as 
l  +  c  or  (  +  0.  An  alternative  way  of  determining  the  boundary 
condition  on  the  near  wake  is  to  use  the  velocity  instead  of 
the  pressure, 

V  -r  +  u2  (1-  -§->  •  <3-2' 
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Using  the  same  potential  function  and  cascade  mapping  function 
as  before  (Furuya,  1978) ,  the  physical  flow  field  will  be  mapped 
into  the  potential  plane  (Figure  3.3(b))  and  then  into  the  trans¬ 
form  plane  5  (Figure  3.3(c)).  The  only  difference  in  the  ;-plane 
from  that  of  Furuya  (1978)  is  that  the  new  ^-coordinate  points 
g  and  h  have  been  added,  representing  the  lower  and  upper  near¬ 
wake  end  points,  respectively. 

The  boundary  conditions  for  the  present  new  flow  model  will  be 
expressed  in  terms  of  either  the  real  or  imaginary  part  of  to 
as  follows. 

i)  On  the  Far  Wake: 

x  =  0  ,  -°°<C<g  and  h<£<°°, 

ii)  On  the  Near  Wake: 

£  I 

— )>  for  the  Pressure  Condition, 

C  f 

(3.3) 

l)1* 

— >  for  the  Velocity  Condition, 

C'  (3.4) 


t 
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where  the  Bernoulli  equation  has  been  used  for  the  velocity 
condition  on  the  near  wake.  Since  there  is  no  change  in  the 
branch  cuts  or  singularities,  the  homogeneous  solution  H(;) 
to  the  problem  should  remain  the  same  as  before  (Furuya  1978), 


h(!)  >  yu±i>  . 

The  general  solution  which  satisfies  the  above  mixed-type 
boundary  condition  can  be  readily  written  by  adding  two  new 
terms  to  the  solution  u)(s)  of  Equation  (2.5)  in  Section  2. 
Defining  uw(Z)  as  the  general  solution  to  the  double  wake 
model  flow,  one  can  express  it  in  a  form 


I»U) 


+ 


/u+1)  (S-b)  (c-f) 

V  C-c 


1 

2iTi 


2it 

_ w _  _c 

[HEn  (£'~b>  c 

C'-c 


2it 

_ w _ 

1  +1 ) (£'-b) (S'-f) 
C'-c 


(3.5  ) 


where  <jo(s)  is  given  in  Equation  (2.5)  and 

t  =  Equation  (3.3)  or  (3.4). 
w 


In  this  problem  we  have  a  total  of  9  unknown  quantities,  U2, 
a2.  A,  b,  c,  f,  g,  h,  requiring  9  equations  to  determine  them 
uniquely; 


i)  At  the  upstream  infinity, 
U)  ( —  a^+i  In  y 


ii)  At  the  downstream  infinity, 
<1>(°°)  =  o2 


(2  equations) 


(1  equation) 


iii)  Length  of  first  arc  i 


(1  equation) 


iv)  End  point  of  the  streamline  (C  *  f) 

matches  the  physical  trailing  edge:  (1  equation) 


v)  Continuity  equation 

vi)  Length  of  cavity  matches  the 
specified  length 

vii)  Cavity  end  to  land  on  the  upper 
surface  of  blade 


viii)  <t>.  -  ♦  =  r 

h  g 


(1  equation) 


(1  equation) 


(1  equation) 


(1  equation) 


Application  of  these  conditions  to  Equation  (3.5)  yields  the 
following  nine  equations,  f^w,  i  =  1,  9  where  f  ^ 
g2  and  gc  are  those  used  in  Furuya  (1978)  or  found  in  Section  2 
The  new  equations,  f^w  are  written  by  simply  adding  the 
new  terms  due  to  the  change  for  the  boundary  conditions  to 
these  old  terms. 


i)  At  Upstream  Infinity: 


U,+l)  U,-b)  c,-f)  t 


f) 


(5+D  (b“5)  (C-f) 


C-c 


+  / - £ - 

/y/rr+iTTB^-nT^ 


r^rj +  92(e) 


(3.13) 


In  order  to  determine  the  second  arc  length  S2»  the  end  point 
of  cavity  should  be  known.  The  cavity  shape  is  expressed  in 
terms  of  ^-coordinate; 


x  -  x. 


V  1+0 


=  /c°s  9ot(5')  «• 


(3.14) 


y  -  Yt 


/  sin  W5'1  ar-  «' 


(3.15) 


where 


(5+1) (S-b) (f-5) 
c-5 


H/ 


'+!)  (5'-b)  (T-f) 
5’-c 


dll 

5'-5 


f(5'+l)  (S'-b)  (5 *-f ) 
T^c 


(3.16) 


xBr  yB  are  the  Physical  coordinates  of  the  blade  leading  edge. 
The  end  point  of  the  cavity  (x  ,  y  )  is  obtained  by  setting 
the  upper  bound  of  the  above  integrals  (3.14)  and  (3.15)  to 
be  c.  Once  (x  ,  y  )  are  determined,  the  second  arc  length  S2 
is  calculated  based  on  the  upper  body  profile. 


Equations  (3.6),  (3.7),  (3.8),  (3.9),  (3.11),  (3.17),  (3.18), 

(3.19)  and  (3.20)  now  provide  nine  independent  relationships 
for  nine  unknown  parameters. 

The  same  numerical  method  as  before  (Furuya  1978)  will  be  applied 
to  the  above  nine  functional  nonlinear  equations  for  determining 
the  nine  solution  parameters . 
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3.2 


NUMERICAL  RESULTS 


A  substantial  effort  has  been  taken  for  modifying  the  previous 
computer  program  (PCASLE)  in  order  to  implement  the  theory 
described  above.  The  new  computer  program  designated  "PCASLW" 
has  the  nine  unknown  parameters  for  the  nine  independent  rela¬ 
tionships.  The  convergence  of  Newton's  iterative  method  was 
very  slow,  costing  about  $50  of  Lawrence  Berkeley  Laboratory's 
computer  time  to  obtain  results  for  one  data  point.  Due  to  the 
time  and  budgetary  restraints,  only  a  few  representative  com¬ 
putations  were  made  with  the  new  computer  code  PCASLW. 

Figures  3.4  and  3.5  show  the  lift  and  drag  coefficients  with 
=  4°  for  the  blade  sectional  profile  of  Hydronautics  7607.02 
propeller  at  x  =  .5,  in  comparison  with  those  of  the  previous 
computer  program  (PCASLE) .  The  relationship  between  the  length 
of  cavity  and  cavitation  number  for  the  same  case  is  shown  in 
Figure  3.6.  It  is  seen  from  these  figures  that  the  new  theory 
has  provided  a  significant  improvement  for  both  force  coef¬ 
ficients  and  cavity  length-cavitation  number  relationships. 

A  smooth  connection  between  the  partially  cavitating  regime 
and  supercavitating  regime  was  provided  by  following  the 
result  of  the  new  p/c  cascade  theory  up  to  o  =  .22,  for  this 
case,  and  then  going  into  the  supercavitating  range.  In  doing 
this  the  data  for  a  larger  than  .22,  calculated  with  the  s/c 
cascade  theory,  should  be  discarded  for  the  previously  men¬ 
tioned  reason  (see  Section  2.2). 

However,  the  present  PCASLW  experienced  difficulty  in  obtain¬ 
ing  the  convergent  solution  for  the  case  of  a1  =  8"  at  which 
the  cavity  boundaries  cleared  the  blade  trailing  edge  for  all 
a ' s  at  the  s/c  range  and  thus  a  rather  smooth  matching  was 
obtained  with  the  previous  p/c  cascade  theory  (see  Figures 
2.4  to  2.6).  The  same  difficulty  existed  for  the  blade 
sectional  profile  at  the  spanwise  location  of  60%  of  Hydro¬ 
nautics  7607.02  at  ■  4°.  It  might  be  coincidental  to 


i  mmmmmm  ••■■vs 


observe  that  for  these  cases  the  cavity  boundaries  clear  the 
blade  trailing  edge  for  all  a ' s  in  the  supercavitating  range 
and  the  discrepancy  for  CL  and  CD  at  *-c/c  =  1  is  minor  (see 
Figures  2.7  to  2.9). 

For  the  case  of  *  2"  at  x  =  .6,  although  convergent  solu¬ 
tions  were  obtained  with  the  new  method,  the  calculated  data 
showed  a  peculiar  behavior  similar  to  that  occurring  in  =  2° 
at  x  *  .5;  the  calculated  cavitation  number  decreased  as  the 
length  of  cavity  decreased.  This  type  of  behavior  in  cavity 
length  vs.  cavitation  number  is  anomalous  in  the  cavity  flow 
unless  this  theoretical  finding  is  interpreted  as  an  indication 
for  unsteady  flow  phenomena.  As  a  matter  of  fact.  Figure  2.11 
indicated  that  the  upper  cavity  boundary  even  for  a  relatively 
long  cavity,  i.e.,  i,c/c  =  1.4  chords,  barely  cleared  the  upper 
blade  surface.  It  is  thus  questionable  if  the  partial  cavity 
can  exist  at  such  a  small  incidence  angle  =  2°. 

It  must  also  be  mentioned  that  this  type  of  reversal  relation¬ 
ship  for  the  cavity  length  to  cavitation  number  took  place  for 
the  partially  cavitating  cascade  flow  around  the  plano-convex 
blade  (see  Furuya  1978) . 
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4. 


SUMMARY  AND  RECOMMENDATION 


The  objective  of  the  present  research  work  is  to  incorporate 
the  partially  cavitating  cascade  data  into  the  performance 
prediction  method  for  cavitating  propellers  and  to  determine 
if  any  improvement  will  be  made  on  the  previously  obtained 
results  (Furuya  1978b)  .  The  partially  cavitating  cascade 
theory  was  developed  in  the  first  phase  of  the  work,  the  com¬ 
puter  codes  developed  then  being  called  PCAS  and  PCASL  (Furuya 
1978) .  It  was  found  that  this  theory  can  be  used  only  for 
cascade  blades  having  small  camber  such  as  flat  or  plano¬ 
convex  blades.  For  practical  blade  profiles  having  rather 
high  camber  such  as  those  used  for  supercavitating  propellers, 
it  was  discovered  that  the  calculated  upper  cavity  boundary 
crossed  the  blade  surface,  which  is  physically  unacceptable. 

In  order  to  correct  this  deficiency,  an  additional  boundary 
condition  that  the  upper  cavity  end  point  be  placed  on  the 
blade  surface  was  imposed  on  the  previous  theory  but  by  allow¬ 
ing  the  wake  to  open  at  downstream  infinity.  The  computer 
codes  developed  based  on  this  flow  model  were  designated  as 
PCASE  and  PCASLE.  Although  the  problem  of  the  upper  cavity 
boundary  crossing  the  blade  surface  was  resolved,  a  new  problem 
of  mismatching,  as  the  length  of  cavity  approaches  that  of 
chord,  has  emerged.  It  was  considered  that  one  of  the  main 
reasons  for  mismatching  was  attributable  to  the  difference  in 
the  base  pressure  between  the  partially  cavitating  and  super¬ 
cavitating  cascade  theories,  leading  to  the  development  of 
another  new  theory.  The  new  flow  model  used  here  was  desig¬ 
nated  as  "double  wake"  flow  model  in  which  the  near  wake  was 
formed  immediately  after  the  blunt  trailing  edge  and  followed 
by  the  far  wake.  The  near  wake  pressure  p  was  determined  by  a 
linear  interpolation  between  the  cavity  pressure  pc  and  far 
wake  pressure  p2  as  a  function  of  the  cavity  length.  By  using 
this  interpolation  method  it  was  ensured  that  the  base  pres¬ 
sure  approaches  to  the  cavity  pressure  as  the  length  of 
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cavity  becomes  that  of  chord.  The  new  computer  program,  now 
the  third  one,  called  "PCASLW"  was  applied  to  the  same  blade 
sectional  profiles  of  Hydronautics  7607.02  propeller.  For  a 
certain  combination  of  the  flow  configuration,  the  calculated 
data  smoothly  connected  with  those  of  the  supercavitating 
regime  as  a  function  of  cavitation  number.  It  was  found,  how¬ 
ever,  that  the  program  "PCASLW"  had  difficulty  in  convergence 
for  large  incidence  angle  cases.  This  was  not  surprising 
because  for  such  cases  a  sudden  jump  of  the  cavity  closure 
point  in  the  y-direction  existed  between  the  supercavitating 
and  partially  cavitating  flows. 

Another  problem  with  "PCASLW"  was  that  the  relationship  between 
the  cavity  length  and  cavitation  number  was  reversed  to  that 
usually  found  in  the  cavity  flow,  i.e.,  the  calculated  cavita¬ 
tion  number  decreased  as  the  specified  cavity  length  decreased. 
This  may  be  interpreted  as  an  indication  for  a  possible  unsteady 
flow  phenomenon  or  for  nonexistence  of  partially  cavitating 
flows  for  such  small  incidence  angles. 

The  lift  and  drag  forces  calculated  based  on  PCASE  and  PCASLE 
programs  were  incorporated  into  the  propeller  program  "SCPR0P2" 
for  prediction  of  the  high  speed  propeller  performance  operating 
at  the  partially  cavitating  regime.  Although  the  data  were 
known  to  be  somewhat  inaccurate,  the  intention  of  the  study 
was  to  seek  a  possibility  for  improving  the  prediction  capa¬ 
bility  of  the  SCPROP2  program  for  the  p/c  regime.  In  fact, 
the  results  of  computations  revealed  such  an  improvement 
despite  the  inaccurate  data  employed. 

The  following  conclusions  and  recommendations  based  on  the 
findings  from  the  current  study  will  be  thus  given. 

i)  The  partially  cavitating  cascade  flow  is  much  more  com¬ 
plex  than  its  counterpart,  i.e.,  supercavitating  cascade 
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flow,  due  to  the  fact  that  the  upper  surface  of  the  blade 
comes  into  a  part  of  the  boundary  conditions.  Three  dif¬ 
ferent  theories  and  their  computer  programs  have  been 
developed  to  date  under  the  GHR  project,  but  no  single 
theory  can  completely  cover  various  flow  configurations 
and  cascade  geometries.  Selection  of  an  appropriate 
one  among  these  three  methods  requires  careful  examina¬ 
tion  of  not  only  the  flow  configuration  to  be  solved  but 
also  that  of  the  supercavitating  condition.  No  general 
selection  guideline  exists  except  for  the  fact  that  the 
single  wake  model  (PCASE  and  PCASLE)  provided  better 
results  in  convergence  as  well  as  the  calculated  data 
for  the  cases  in  which  the  cavity  boundary  calculated 
from  the  s/c  theory  clears  the  upper  trailing  edge  near 
£c/c  =  1.  For  the  cases  in  which  the  cavity  boundary 
crosses  the  upper  trailing  edge  near  ^c/c  =  1,  the  double 
wake  model  (PCASLW)  was  superior. 

ii)  Because  of  the  imperfect  state  of  the  three  partially 

cavitating  cascade  theories  developed  here,  it  is  recom¬ 
mended  that  a  more  general  theory  should  be  developed 
with  a  different  set  of  boundary  conditions. 

iii)  In  order  to  develop  a  better  flow  model,  more  detailed 
information  is  necessary  from  experiments.  This  should 
include  the  relationship  between  the  cavity  length  and 
cavitation  number,  forces,  base  pressure  and  stability 
of  the  cavity  flow  for  various  flow  and  blade  configura¬ 
tions. 

iv)  The  accuracy  of  the  three  theories  developed  here  should 
be  determined  only  in  comparison  with  experimental  data. 


v)  Based  on  the  trial  test  conducted  during  this  phase  of 

work,  it  was  indicative  that  more  accurate  partially  cavi- 
tating  cascade  data,  if  available  and  used  in  the  high 
speed  propeller  theory,  will  improve  the  prediction 
accuracy  for  cavitating  propeller  performance.  It  is, 
however,  too  early  at  this  stage  to  determine  which  flow 
model  or  theory  should  be  used  for  this  purpose.  Again 
we  should  await  the  results  of  experiments,  which  will 
clarify  many  uncertainties  mentioned  above. 
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TABLE  2.1 


Lift  and  Drag  Data  Near  A  =  1  Calculated  from  the  Supercavitating 
(s/c)  Cascade  Theory  and  Partially  Cavitating  (p/c)  Cascade  Theory 


i)  At  spanwise  location  x  =  0.5  of  Hydronautics  7609.02 
Propeller 


Theory 

Item  ' 

s/c  Cascade 
Theory  (1975) 

p/c  Cascade 
Theory  (Present) 

Cavitation  Number 

0 

.450 

.495 

p 

II 

00 

0 

Cavity  Length 
l/o 

1.020 

.900 

CL 

.866 

.871 

CD 

.181 

.138 

j 

a 

.220 

.201 

*'3 

l/c 

1.091 

.950 

CL 

.405 

.369 

■ 

CD 

.0546 

.0300 

TABLE  2.1  (Con't.) 


x  =  .6 


^"^-^Theory 

Item 

s/c  Cascade 
Theory 

p/c  Cascade 
Theory 

0 

.200 

.284 

l/c 

1.190 

.900 

CL 

.398 

.505 

CD 

.0954 

.0421 

0 

.115 

.102 

l/c 

1.400 

.900 

CL 

.261 

.231 

CD 

.0246 

.0126 

_ 


L 


* 


t 
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Figure  2.1  Cavity  Boundary  Calculated  with  the  Previous  P/C  Cascade 
Theory  (Furuya  1978)  for  Hydronautics  7607.02  Propeller 
Blade  Section  at  50%  Span  Position,  a.,  =  4°  and  i  =  . 35 


LIFT  COEFFICIENT 


L  10-1  I _ _ 1 _ 1 - U  1  1  OO.lJ 


LIFT  COEFFICIENT 


DRAG  COEFFICIENT 


Figure  2.8  CD  vs.  o  for  the  Same  Blade  Section 
as  that  of  Figure  2.7 


Figure  2.10  Locations  of  the  Calculated  Cavity  Boundaries  for 
Blade  Profile  at  the  Spanwise  Position  of  50%  of  Hydronautics  7607.02 


Figure  2.11  The  Same  as  that  of  Figure  2.10 
except  that  the  Spanwise  Location  is  60% 


Figure  2.12  Comparison  of  the  Present  Results  for  K^,  and  n  of 

Hydronautics  7607.02  at  oVa  =  .343  with  those  of  the 

Previous  Theory  (Furuya  1978b)  and  Experiments  (Bohn  1977 
and  Peck  1977) 


(a)  Partially  Cavitating  Condition 


(b)  Supercavitating  Condition 


Figure  3.1  Comparison  of  the  Flow  Configuration 
Between  the  Partially  Cavitating  Con¬ 
dition  and  Supercavitating  Condition 
as  the  Cavity  Length  approaches  to  the 
Chord  Length 
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closed  at  infinity 


3  Schematic  Diagram  for  the  Double 
Wake  Model  for  the  Partially 
Cavitating  Cascade  Flow  and  its 
Mapping  Planes 


0  .1  .2  .3  .4  .5  .3  .7 

CAVITATION  NO.  t 


Figure  3.4  Lift  Coefficient  vs.  Cavitation  Number 
Calculated  with  the  New  Partially  Cavitating  Cascade 
Theory  in  Comparison  with  the  Previous  Theory  at 
x  *  0.5  for  Hydronautics  Propeller  7607.02 


tun  on  s/e  an euc  ,n«tc*r 


Figure  3.5  Drag  Coefficient  vs.  Cavitation 
Number  for  the  Same  Case  as  Figure  3.4 


to 

6 


to 

6 


o 


54 


Figure  3.6  Cavity  Length  vs.  Cavitation  Number  for 
the  Same  Case  as  Figures  3.5  and  3.6 


